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In contrary to other ID momentum-conserving lattices such as the Fermi-Pasta-Ulam /3 (FPU-/3) 
lattice, the ID coupled rotator lattice is a notable exception which conserves total momentum while 
exhibits normal heat conduction behavior. The temperature behavior of the thermal conductivities 
of ID coupled rotator lattice had been studied in previous works trying to reveal the underlying 
physical mechanism for normal heat conduction. However, two different temperature behaviors 
of thermal conductivities have been claimed for the same coupled rotator lattice. These different 
temperature behaviors also intrigue the debate whether there is a phase transition of thermal con¬ 
ductivities as the function of temperature. In this work, we will revisit the temperature dependent 
thermal conductivities for the ID coupled rotator lattice. We find that the temperature dependence 
follows a power law behavior which is different with the previously found temperature behaviors. 

Our results also support the claim that there is no phase transition for ID coupled rotator lattice. 

We also give some discussion about the similarity of diffusion behaviors between the ID coupled 
rotator lattice and the single kicked rotator also called the Chirikov standard map. 

PACS numbers: 05.60.-k,44.10.-|-i,05.45.-a 


I. INTRODUCTION 

The exploration of underlying mechanism for anoma¬ 
lous and normal heat conduction in low dimensional sys¬ 
tems represents a huge challenge in the area of statisti¬ 
cal physics After enormous efforts for more than 

one decade from numerical simulations [sl-ld^l . theoretical 
predictions [45l - [57l| and experimental observations [ssl - l^ . 
there is still no consensus for the exact physical mecha¬ 
nism causing anomalous heat conduction. It is believed 
that momentum conservation plays an important role 
in determining the actual heat conduction behavior. In 
general, ID non-integrable lattices with momentum con¬ 
serving property should have anomalous heat conduction 
where the thermal conductivity k diverges with the lat¬ 
tice size TV as K oc N°’ where 0 < a < I M. However, 
the ID coupled rotator lattice is a well known exception. 
It exhibits normal heat conduction behavior despite its 
momentum conserving nature [6l|, . 

The normal heat conduction in ID coupled rotator 
lattice was discovered via numerical simulations by two 
groups independently [^ . In order to understand 
the underlying mechanism, the temperature dependence 
of thermal conductivity has been studied in detail in 
both works. In Ref. [^, the temperature dependence of 
thermal conductivity of coupled rotator lattice has been 
found to be like k{T) oc where AU is a posi¬ 

tive constant. However, in Ref. [b^, a temperature de- 
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pendence of thermal conductivity of k{T) oc was 

claimed where A is a positive constant. There it was 
also argued that a possible phase transition at temper¬ 
ature around T ^ 0.2 — 0.3 where heat conduction is 
normal above this temperature while anomalous below 
this temperature exists [b^ • Although this phase transi¬ 
tion statement was challenged as a finite size effect [b^, 
a later work supported this phase transition conjecture 
after deriving a similar temperature dependent thermal 
conductivity as k{T) oc (5flj |. 

Most recently, the ID coupled rotator lattice re¬ 
attracts much attention due to the new finding of simul¬ 
taneously existing normal diffusion of momentum as well 
as heat energy (b4| . It is argued that the normal be¬ 
havior might be due to the reduced number of conserved 
quantities which is unique for periodic interaction poten¬ 
tials [110. But the new debate is that whether the 
stretch is conserved or not in ID coupled rotator lattice. 
The investigation of ID coupled rotator lattice will be 
the key to unravel the true mechanism behind the con¬ 
nection between momentum conservation and normal or 
anomalous heat conduction. Therefore, it is the right 
time to revisit the temperature dependence of thermal 
conductivity of ID coupled rotator lattice as a first step. 

In this work, we will revisit the temperature depen¬ 
dence of thermal conductivities for the ID coupled rota¬ 
tor lattice. We find that the temperature dependence is 
neither k(T) oc nor K(r) oc as previously 

claimed [bl|, |b2| . The actual temperature dependence is 
a power-law dependence as k(T) oc The possible 

connection with the momentum diffusion of single kicked 
rotator or the Chirikov standard map has also been dis¬ 
cussed. In order to determine whether there is a phase 
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transition, we also present the temperature dependent 
thermal conductance at different system sizes. All the 
thermal conductances for different sizes collapse to the 
same value at low temperatures while approach to the 
power-law behavior as k(T) oc T~^-^ at high tempera¬ 
tures. However, the crossover temperature decreases as 
the system size increases. This fact indicates that there is 
no phase transition between normal and anomalous heat 
conduction. In thermodynamical limit, the heat conduc¬ 
tion is normal for all temperatures except the trivial zero 
temperature point. In Sec. II the lattice model will be 
introduced. Numerical results and discussions will be 
presented in Sec. Ill and the results will be summarized 
in Sec. IV. 


II. MODEL 


The Hamiltonian for the ID coupled rotator lattice is 
defined as the following [^, : 


N 


^ = E 


+ cos {qi+i - Qi)] 


( 1 ) 



FIG. 1; (color online). Temperature profiles for the ID cou¬ 
pled rotator lattice at (a) T = 0.1 and (b) T = 0.8. The 
lattice sizes are N = 50,100, 200, 400, 800 and 1600 and the 
color rule of the lines are the same for both (a) and (b). The 
first and last atom are put into contact with a Langevin heat 
bath with temperature set as = T(l± A) where A = 0.1 
here. 


where qi and pi denote the displacement from equilibrium 
and momentum for f-th atom, respectively. The mass of 
the atom m and the Boltzmann constant ks has been set 
into unity. The parameter K with energy dimension rep¬ 
resents the coupling strength of the inter-atom potential. 
Therefore, the system temperature T can be rescaled by 
T jK and we can also set AT = 1 for simplicity [b^ ■ The 
equations of motion for i-th atom are 


qi = Pi 

Pi = K sin (gi+i - qi) - K sin {qi - qi^i) (2) 


At low temperature limit, the displacements are small 
values so that \qi+i — qi\ 1. The Hamiltonian can be 
expanded into the Taylor series as: 


H 


E 


^ - Qi)'^ 

2 ^ 2 


{q^+l - qi)'^ {q,+i - qif 


4! 


b! 


( 3 ) 


This Hamiltonian will approach to the integrable Har¬ 
monic lattice only at zero temperature T = 0. The first 
stable nonlinear potential term will be the sextic poten¬ 
tial. 

On the other hand, at infinite high temperature limit 
T = oo, the rotator lattice will approach to another inte¬ 
grable system consisting of N independent and isolated 
free particles as 


N 


^ 2 


( 4 ) 


This is because the kinetic energy is proportional to the 
temperature as (pf) = T and the potential energy in Eq. 
o is confined by the cosine function. 

Without getting into numerics, we can get a qualitative 
picture of the thermal conductivity of ID coupled rota¬ 
tor lattice as the function of temperature. The thermal 
conductivity k{T) will diverge in the low temperature 
limit approaching to the harmonic limit and will decay 
to zero in the high temperature limit as approaching to 
the unconnected N free particles. As a general picture, 
the thermal conductivity k{T) will decrease as the tem¬ 
perature increases. 

In non-equilibrium numerical simulations, we put the 
first and last atom into contact with Langevin heat bath. 
The temperatures for left and right heat bath are set as 
Tl/r = T{1 ± A) where T denotes the average tempera¬ 
ture and A = O.I gives rise to the temperature gradient 
along the lattice. The fixed boundary conditions with 
qo = qN+i = 0 are also been applied. 


III. RESULTS AND DISCUSSIONS 

Before we discuss the results of thermal conductivi¬ 
ties, we first show the temperature profiles. In Fig. [TJ 
the temperature profiles at two different temperatures 
T = O.I and T = 0.8 are plotted. The lattice sizes are 
chosen a,s N = 50,100,200,400,800 and IbOO. For rela¬ 
tive low temperature at T = O.I, the temperature jumps 
at two boundaries are obvious. However, the tempera¬ 
ture jumps are reduced with the increase of lattice size 
N as can be seen in Fig. [IJa). In Fig. [T](b) where the 
temperature is relatively high at T = 0.8, all the tern- 
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FIG. 2: (color online). Thermal conductivities n{T) as 
the function of temperature for different lattice sizes N = 
50,100, 200, 400, 800 and 1600. The straight line is propor¬ 
tional to T~^'^ which describes the temperature behavior of 
K(r) at high temperature region very well. 


FIG. 3: (color online). Thermal conductivities fi^(T) as the 
function of temperature for lattice size N = 1600. The data 
are taken from Fig. [2]and can be best fitted to be a power-law 
dependence as k{T) oc 


perature profiles collapse to the same straight line as the 
temperature jumps are very small for all lattice sizes. 

In order to obtain the temperature dependence of ther¬ 
mal conductivities, we need first define the way how k{T) 
can be calculated numerically. We notice that the tem¬ 
perature profiles all collapse to the same straight line if 
the temperature jumps can be ignored at high tempera¬ 
tures or long lattice sizes. It is appropriate to dehne the 
thermal conductivity k{T) as: 


<T) 


JN 

Tl-Tr 


( 5 ) 


where J = (Ji) is the average heat flux along the lattice in 
the stationary state and the local heat flux Ji is defined 
as Ji = —piK sin {qi^i — qi) via the energy continuity 
equation. 

In Fig. [21 the thermal conductivities k{T) as the func¬ 
tion of temperature are plotted for different lattice sizes 
N = 50,100, 200,400,800 and 1600. At high temperature 
region, all the thermal conductivities k{T) for different 
lattice sizes collapse together indicating the saturation 
of thermal conductivities as the increase of lattice size. 
This is characteristic for lattices with normal heat con¬ 
duction. As the temperature decreases, the thermal con¬ 
ductivities first increases and then becomes flat. This is 
because the phonon mean free paths are getting longer 
as the temperature reduces and the ballistic regime will 
be approached if the phonon mean free paths are longer 
than the lattice size. The boundary jumps will dominate 
the temperature profiles as in Fig. (Ha) and the defini¬ 
tion of k{T) of Eq. ([5]) will no longer reflect the actual 
thermal conductivities. 

At high temperatures, it is clearly seen that the ther¬ 
mal conductivities follows a power-law dependence as 
k(T) oc For the longest size we considered here as 

N = 1600, this power-law behavior can be best fitted for 


more than two orders of magnitudes for the k{T) value 
as shown in Fig. O This also explains the poor fitting 
of the k{T) oc dependence in Ref. and the 

K(r) oc dependence in Ref. 

This power-law dependence of k{T) oc T~^ ‘^ cannot be 
explained by the effective phonon theory which is able 
to predict the temperature dependent thermal conduc¬ 
tivities for other typical ID nonlinear lattices such as 
FPU-/3 lattice and generalized nonlinear Klein-Gordon 
lattices [bSl - f^ . According to the effective phonon the¬ 
ory, the thermal conductivity for low temperature rotator 
lattice with Hamiltonian of Eq. m can be derived as 

k{T) oc i oc T-2 (6) 

e 

where e is the nonlinearity strength with the following 
temperature dependence 


((g»-n - 

{(g*+i - qi)^) 


ocT^ 


( 7 ) 


at low temperature region. Here we consider the sextic 
potential term as the lowest nonlinear term because the 
dynamics governed by the negative quartic potential term 
is unstable. Therefore, the actual temperature behavior 
of k{T) oc is steeper than the prediction of k{T) oc 

T~^ from effective phonon theory. 

Although it is difficult to unravel the exact physical 
mechanism behind the power-law dependence of ther¬ 
mal conductivities for coupled rotator lattice, it is very 
helpful to look into the transport properties of the single 
kicked rotator (the Chirikov standard map) [t^. As the 
name indicates, the coupled rotator lattice is a kind of 
connected kicked rotators. The equations of motion for 
the single kicked rotator are 


Pn+I-Pn = Ksm{qn) 
Qn-\-l Qn — Pn+1 


( 8 ) 
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FIG. 4; (color online). Momentum diffusion constant Dp as 
the function of temperature for ID coupled rotator lattice. 
The numerical data are obtained via equilibrium MD simula¬ 
tions as in Ref. [^. The solid line of T~^'^ is guided for the 
eye. 

where and pn denote the coordinate and momentum 
after n-th kick. 

The variation of momentum p is unbounded and can 
be characterized by a normal diffusion [tI, [t^] 

(Ap2(t)) - Dt (9) 

where D is the diffusion constant. This is similar to the 
normal momentum diffusion for ID coupled rotator lat¬ 
tice. According to Ref. [zilzS , the diffusion constant D 
depends on the coupling strength K as 

D oc {K-1.2f, 1.2 <K <4 (10) 

D (X K^, K >4 (11) 

where 1.2 is the chaos threshold for the kicked rotator. 

If one assume the transport properties of coupled rota¬ 
tor lattice are the same as that of single kicked rotator, 
one would expect that the energy diffusion constant De 
as well as the momentum diffusion constant Dp oi cou¬ 
pled rotator lattice should also follow the same depen¬ 
dence as that of single kicked rotator as in Eq. (flUl) . To 
translate the K dependence into T dependence, we notice 
that the parameter K in single kicked rotator of Eq. (|S]) 
plays the same role as that in the coupled rotator lattice 
as in Eq. ©■ As we have discussed above, the dynamics 
of coupled rotator lattice can be rescaled with T/K. The 
K dependence should be inversely proportional to the T 
dependence. Low K value region corresponds to the high 
T region and vice verse. And for normal heat conduction, 
the thermal conductivity n is proportional to the energy 
diffusion constant Dp- This will finally give rise to the 
prediction of temperature dependent thermal conductiv¬ 
ity k(T) of coupled rotator lattice from the analogy with 
single kicked rotator: 

k(T) oc T-^, lowT (12) 

k(T) oc T-3, highT (13) 


FIG. 5: (color online). Heat conductance cr as the 
function of temperature T for different lattice size N = 
50,100, 200,400, 800 and 1600. All the other parameters are 
the same as in Fig. [2] The left-pointing arrow represents the 
trend of decreasing crossover temperature as the lattice size 
increases. 


As a consistence check, we plot the momentum diffu¬ 
sion constant Dp for the coupled rotator lattice as the 
function of temperature in Fig. |4l The same temperature 
dependence of Dp oc has been obtained. 

Therefore, the prediction of k(T) oc T~^ with anal¬ 
ogy to the single kicked rotator at high temperature or 
low K region is close to the numerical observation of 
k{T) oc This indicates there might be some deeper 

connection between the diffusion behavior of single kicked 
rotator and ID coupled rotator lattice. In addition, the 
analogy analysis also predicts a k(T) oc T~^ behavior at 
low temperature region which is the same as the predic¬ 
tion from the effective phonon theory of Eq. (|6l). The rea¬ 
son why this temperature behavior cannot be observed in 
numerical simulations might be due to the severe finite 
size effect at low temperature region as can be seen in 
Fig. [Ha). 

In the final part we will discuss the issue about the 
possible phase transition between normal and anomalous 
heat conduction for ID coupled rotator lattice. As we 
have shown in Fig. [31 the actual temperature dependence 
of k{T) is a power-law dependence as k{T) oc This 

indicates the analytic derivation of k{T) oc in Ref. 

[ 5 ^ originally for a lattice with pinned on-site potential 
can not be extended to ID coupled rotator lattice. 

On the other hand, the claim that the transition tem¬ 
perature is around (0.2 — 0.3) in Ref. has been chal¬ 
lenged in Ref. with numerical simulations of longer 
sizes. The possible phase transition could be a finite size 
effect. Here we reconfirm the finite size effect by giving 
a more illustrative picture in Fig. [5] The heat conduc¬ 
tance a = k(T)/N as the function of temperature has 
been plotted for different lattice sizes. At high temper¬ 
ature regions, all the heat conductances follow the same 
temperature behavior of cr oc As the tempera- 
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ture decreases, the phonon mean free path will overcome 
the lattice size and the heat conductance will saturate 
to a value determined by the harmonic limit of coupled 
rotator lattice. As can be seen from Fig. [SJ the heat 
conductance for short lattice will first bend and become 
flat as the temperature decreases. This crossover tem¬ 
perature decreases as the lattice size increases smoothly 
and no sign of phase transition can be observed. It is 
also noticed that the crossover temperature happens to 
be around (0.2 — 0.3) for lattice sizes of a few thousands. 

IV. SUMMARY 

We have systematically investigated the temperature 
dependence of thermal conductivities of ID coupled ro¬ 
tator lattice. The actual temperature dependence is a 
power-law dependence of k{T) cx T~^'^ which is different 
with the observations of previous studies. The possible 


connection with the single kicked rotator or the Chirikov 
standard map has been discussed where a k,{T) oc T~^ 
dependence can be implied. Our results also recon¬ 
firm that the previously claimed possible phase transition 
should be a finite size effect. 
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